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Abstract. We construct a generalization of the multiplicative prod- 
uct of distributions presented by L. Hormander in [L. Hormander, The 
analysis of linear partial differential operators I (Springer- Verlag, 1983)]. 
The new product is defined in the vector space _4(R) of piecewise smooth 
functions / : R — *■ C and all their (distributional) derivatives. It is as- 
sociative, satisfies the Leibniz rule and reproduces the usual pointwise 
product of functions for regular distributions in A(M). Endowed with 
this product, the space _4(R) becomes a differential associative algebra 
of generalized functions. By working in the new ^4(M)-setting we deter- 
mine a method for transforming an ordinary linear differential equation 
with general solution t[> into another, ordinary linear differential equa- 
tion, with general solution xn^i where xn is the characteristic function 
of some prescribed interval fi C R. 

1. Introduction 

Two of the most interesting properties of the space of Schwartz distri- 
butions Z>'(R n ) over M. n are that the space of continuous functions C°(M n ) 
is a subspace of T>'(M n ) and that differentiation is an internal operation in 
T>'(R n ) [221 [261 H]. On the other hand, its major limitation is that it only 
displays the structure of a vector space and not that of an algebra. This has 
been known since 1954 when L. Schwartz proved that there is no associative 
and commutative algebra of generalized functions (^(K"), + , o) satisfying 
the three following properties: 

(1) The space of distributions D'(R n ) over M. n is linearly embedded into 
_4(IR n ) and f(x) = 1 is the identity in 



(2) The restriction of the product o to the set of continuous functions C°(IR n ) 
reproduces the pointwise product of functions. 

(3) There exist linear derivative operators -J^r : *4(M n ) — ► „4(M n ), (i = 
1, .., n) that: (a) Satisfy the Leibniz rule and (b) Their restrictions to V^W 1 ) 
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coincide with the usual distributional derivatives. 

This theorem became known as the Schwartz impossibility result [21]. In 
spite of it, certain products of distributions emerge rather naturally in sev- 
eral fields of research in both mathematics and physics. Common examples 
can be found in quantum electrodynamics and particle physics when deal- 
ing with self interacting terms; in hydrodynamics in the formulation of the 
dynamics of shock waves; and, more generally, whenever one wants to con- 
sider differential equations where either the coefficients or the prospective 
solutions are non-smooth. 

These issues motivated several proposals towards a precise construction 
of a multiplicative product of distributions. Among these, the most famous 
example is due to J. Colombeau [U[2j[3], who constructed differential alge- 
bras of generalized functions Q(M n ) satisfying conditions (1) and (3) with 

(2) holding only in C°°(]R r! ') (another proposal of a globally defined product 
of distributions can be found in |19t I20j). Colombeau algebras are a fine so- 
lution for the problem of constructing associative and commutative algebras 
G(M. n ) in the situation: 

(1) C°°(R n ) C C7°(R n ) C V'(R n ) C g(M n ) 

and satisfying the maximal number of properties stated in the Schwartz 
impossibility result. Its main disadvantage might be that, in practical ap- 
plications, one is required to leave the simpler distributional framework and 
work in the more involved setting QiW 1 ). This is because, in general, the 
Colombeau product of two distributions is no longer a distribution. More- 
over, C°(R n ) is not a subalgebra of C?(R n ). 

In this paper we want to consider the following problem: Since it is not 
possible to construct a superspace of T>'(M. n ) satisfying the properties (1) to 

(3) stated in the Schwartz impossibility result, one could, at least, try to 
construct (the largest possible) subspace of P'(R n ) that satisfies properties 

(2) and (3). More precisely: 

Problem 1 

Determine associative (possible non-commutative) algebras (,4(R n ), +, *) 
satisfying 

(2) C°°(R n ) C A(M n ) C 2?'(R n ) 

and such that: 

(i) The restriction of the product * to the set A(M n ) n C°(R n ) coincides 
with the usual pointwise product of functions, i.e. / * g = fg, for all 
f,g e A(M n )nC°(R n ). 

(ii) There exist derivative operators in ^4(R n ), which are of the form : 
„4(R n ) -» A(M n ) and: 
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(ii-1) coincide with the restrictions to A(M. n ) of the usual distributional 
derivatives in V'(M n ), 

(ii-2) satisfy the Leibniz rule. 

In this paper we will address this problem in one dimension. We will take 
•A(R) to be the space Cf°(K) of piecewise smooth functions together with 
their (distributional) derivatives to all orders, and define a new product * 
such that (.A(IR), +, *) becomes an associative (but noncommutative) algebra 
satisfying properties (i) and (ii) above. The new product * is a generalization 
of the product of distributions with non-intersecting singular supports that 
was proposed by L. Hormander in [pag. 55 of ref.[12j]. Our task here 
will consist of extending the domain of Hormander's product to include the 
case of distributions with intersecting singular supports. Clearly, this is not 
possible for the entire set P'(R), but we will prove that it is possible in 
A(R). 

It is important to remark that: 

(a) A(M) is not a subalgebra of <7(R) because the product ★ is not the re- 
striction of the Colombeau product o to ,4(R). In fact, Colombeau's product 
is not an inner operation in A(M) and it does not even satisfy / o g £ D'(R) 
for all f,g& A(M.). 

(b) In many specific problems, where products of distributions are present, 
what is at stake are products involving piecewise C 00 (M)-functions and dis- 
tributional derivatives of these functions. This is the typical case when con- 
sidering differential equations with distributional coefficients or non-linear 
terms. For these cases, the algebra ,4(R) displays almost optimal proper- 
ties. It provides a sufficiently general setting without ever leaving the space 
£>'(R). 

(c) The main limitation of A(M) is that it is defined over M., only. Moreover, 
the details of the derivation of the product * indicate that the extension of 
its domain to functionals of several variables might not be easy. We hope 
to study this issue in a forthcoming paper. 

In the second part of this paper, we use our product of distributions to 
address the following problem: 

Problem 2 

Let C 1 be an open interval and let be the characteristic function of 
0, (i.e. xn(x) = 1 for x € f2 and zero otherwise). 

Let us consider a generic, linear and ordinary differential equation and let 
ip € C°°(M) be its general solution on R. 

The problem we want to address is that of deriving a new, also linear and 
ordinary differential equation, but displaying the general solution xoV'- 

The crux of the matter here is that it is the general solution (i.e. irrespec- 
tively of the boundary conditions) that has to be confined to the interval 
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ft. The case where we only impose the confinement of a particular solu- 
tion can be solved by adding to the original differential equation a suitable 
non-homogenous distributional term. Differential equations of this type and 
their relation with the Green's function theory and with functions with jump 
discontinuities have been investigated in [HI [T8J, [23] . In general the dis- 
tributional term produces a jump discontinuity of fixed strength that is able 
to confine a particular, but not all solutions [14j . In several cases this is 
an important limitation. For instance, if we want to construct a (version 
of a given) differential operator whose eigenfunctions are all confined to an 
interval f2 (a relevant problem, for example, for the global formulation of 
quantum systems with boundaries [13]) the distributional non- homogeneous 
term method is useless. 

By working in A(M) we will provide a general method for constructing the 
differential equations, solutions of Problem 2. Their most interesting fea- 
ture is that they display an extra term written in terms of a distributional 
coefficient. This is why they require the use of a product of distributions. 
The new differential equations yield a global formulation for a class of dif- 
ferential problems defined on bounded domains and may find interesting 
applications in those fields of mathematical physics and dynamical systems 
where global analysis is required [13J. An example of such applications is 
the derivation of a globally defined time independent Schrodinger equation 
for systems with boundaries [131 S El E] an d the related issue of quantum 
confinement in several non-local formulations of quantum mechanics, such 
as the deformation and the De Broglie-Bohm formulations [16} El [6] . These 
problems were studied in [U [TJ [5] using an approach related to the results 
that will be presented here. 

This paper is organized as follows: In section 2 we introduce the notation, 
provide a concise review of some relevant results on Schwartz's distributions 
and define the space «4(M). In section 3 we define a product in A(M) and 
study its main properties. In section 4 we address the problem of construct- 
ing globally defined differential equations displaying solutions of confined 
support. 

2. Schwartz distributions and the space A(R) 

In this section we introduce the notation, review some classical results in 
the theory of distributions and define the space A(M) (in the next section we 
will prove that A(M) is an algebra). Let C M be an open interval, C m {Q) 
the vector space of complex functions on f2 which are m-times continuously 
differentiable, and let C™(Q) be the vector space of piecewise C m -functions 
on f2, i.e. of functions that are m-times continuously differentiable, except 
on a finite set, where they and all their derivatives up to order m have 
finite left and right limits. Let also T>(Q) be the vector space of infinitely 
smooth, complex valued functions with support on a compact subset of O. 
The set T>(£1) is usually named the set of test functions. The set of Schwartz 
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distributions T>'(Q) is the topological dual of X>(0). Its elements are linear 
and continuous functionals F acting on T>(£1) [26\: 

(3) F : V(fl) -> C; t -> F(t) =<F,t> 

where -F(t) and < F,t > are two alternative notations for the action of F 
on t. 

In our notation the lowercase roman letters from the middle of the al- 
phabet (f,g,h...) will be used for arbitrary functions, those from the end 
of the alphabet (t,s,r, ...) will be reserved for test functions and the upper 
case roman letters from the middle of the alphabet (F, G, J, K) will be used 
for distributions. H is the Heaviside step function. By -£-^f and we 
denote the nth-order distributional derivative of /. The notation /', /" will 
also be used for the first and second order derivatives. 

A distribution F is regular if it acts as < F, t >= f f(x)t(x)dx for some 
locally integrable function /. If / is continuous we shall make the identi- 
fication F = f because / is determined by F uniquely. Furthermore, the 
restriction of a distribution F : T>(Q) -> C to an open set 3 C O is the 
distribution F E : P(3) -> C; t ->< F s ,t >=< F,t >. Two distributions 
F,G £ D'(Q) are identical on some open set 5 C iff F-, = G-=. The null 
set of F (denoted null F) is the largest open set 3 where Fs = 0. The 
complement of this set is the support of F, denoted supp F. The singular 
support of F is defined as the complement of the largest open set 3 where 
F = f for some / £ C°°(3). The standard notation for this set is sing supp 
F. Another important notion is that of (weak) convergence in P'(Q). A se- 
quence of distributions (i ? n )i<n<oo is said to converge (weakly) to F G V'(Q) 
iff lim^oo < F n , t >=< F,t> for all t € T>(Q). 

The following definition is very important for the sequel: 

Definition 2.1: The space .4.(0) 

The space of special distributions A(£l) on f2 is the space of piecewise smooth 
functions C^(Q) (regarded as distributions) together with their distribu- 
tional derivatives to all orders. 

Remark 2.2. All functions / £ C^°(0) can be associated with a distri- 
bution through the prescription < f,t >= J f(x)t(x)dx. In this sense, 
C^°(0) C P'(O) and so A(Q) is a subset of V(Q). One easily realizes that 
A(£l) is, in fact, a linear subspace of T>'(fl). It also follows from its definition 
that A(£l) is closed under differentiation. In the next section we will define a 
product -k in A(M) and prove that (*4(M), +,*) is an associative, differential 
algebra. 

Remark 2.3. For each / £ C?°(r2) there is a finite set Vf = {x± < X2 < 
... < x N } C O such that / £ C°°(Q\Vf). Let h =}xi,x i+ i[, i = l,..,N- 1; 
Iq =] — oo, xi[flfi and In =]%Ni +oo[nSl. Since / and all its derivatives dis- 
play finite limits at Xi, i = 1, .., N there is (by Whitney's extension theorem 
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[2U HO]) a collection of functions /, G C°°(U), i = 0,..,N such that their 
restrictions to 1% coincide with those of /, i.e: 

(4) f(x)=fi(x) if x£l l ; i = 0,...,N 
Hence, on Q\Vj we have: 

N 

(5) / = 5>i/i 

where Xi( x ) = 1 if zc £ and Xi(aj) = if x ^ Ij is the characteristic 
function of I{. Conversely, if / is of the form (5) then / G Cp°(0). These 
statements are also valid in distributional sense, i.e. / G C^°(f2) (regarded 
as a subset of - Remark 2.2) iff there is a finite set Vj C Vt and a 

collection of functions fi G C°°(r2), i = 0, .., iV = JjV/ such that, in the sense 
of distributions, / can be written in the form (5). 

Remark 2.4. It is clear that C°°(Q) is a subspace of A(tt), but C°(ft) 
is not. For instance /(x) = \f\x\ ^ >4(M). The Dirac delta distribution 8 
and all its derivatives are elements of A{Q). 

Lastly, we review the content of two well-known theorems, which will be 
important in the sequel. Proofs may be found in [26]. The first one concerns 
the (re)construction of a global distribution from a set of local ones: 

Lemma 2.5. Let N G N and let Q.^ C fi, k = 1,..,N be a finite col- 
lection of open sets covering Q. Let G 2?'(f2fc), k = 1, .., N be a collection 
of distributions such that Fi = Fj in $7j H Qj for all i,j = 1,..,N. Then 
there is one and only one F G T>'(£1) such that Fn k = F^ for all k = 1, .., N. 

The second one is a consequence of the completeness of V{Vt) with re- 
spect to the weak convergence: 

Lemma 2.6. Let F e G T)'{Vt) be a one parameter family of distributions 
with e > 0. If lim (E _ + o+ < F € ,t > exists for every t G T>(£1) then: 

(6) F : V(n) — ► C, t — X F, t >= lim < F e ,t > 

is an element o/X>'(f2). 

3. The algebra of special Distributions and multiplicative 

products 

In this section we will provide two alternative characterizations of the 
space A(M) (Theorems 3.1 and 3.2), study Hormander's product of distribu- 
tions with non intersecting singular supports (Definition 3.4 to Theorem 3.7) 
and make the proposal of the new product * (Definition 3.10 and Theorems 
3.13 to 3.16). These results prove that (A(R),+,*) is indeed a differential 
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associative algebra and a solution of Problem 1. To make the presentation 
simpler, we will assume that = R, but our results are still valid for an 
arbitrary open interval Qcl. 

The next theorem follows from Remark 2.3 and provides an important 
characterization of «4(R): 

Theorem 3.1. F G .A(R) iff there is a finite set of real numbers Vf C M 
and two Schwartz distributions f and A^ F \ such that: 

(7) F = f + A^ 

where f G C£°(R) n C°°(R\Vf), and A^ is of the form: 

N 

(8) A(^) = £ A(f) = Z,C k (w)SM(z-v,) 

weV F w£V F k=0 

for some N G No and Ck{w) G C. 

Proof. If F G A(M), then F is a nth-order (distributional) derivative of 
some function g G C£°(R) (cf. Definition 2.1). It follows from Remark 2.3 
that, in the sense of distributions, g can be written as: 

(9) g = go + ^2 H(x - w)g w 

W&Vg 

where H is the Heaviside step function, go,g w £ C°°(R) and V g is a finite 
set of real numbers. We then have: 

(10) F = ^ = ^+EE (p ) H[n ~ p) ( x - 

W£Vg p=o 

where (™) = ( ra ^i p! is the binomial coefficient. It is easy to realize that F 
can be cast in the form (7) by setting: 

(11) / = ^ 

weVg 

n-l 

weVg p=o 

n-l 

= EE^w^'^-") 

W&Vg P = 

and it is clear that i) / G C£°(R) n C°°(R\V g ) and that ii) A^ is of the 
form (8)). 

Conversely, if F is given by eqs.(7,8) then F is the (iV + l)th-order derivative 
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of: 

(12) 9 = 90+ £ £ W^V HiX ~ W){X ~ ^ 

w&v F k=o v '' 

where: 

( f*dx f(x ) , iV = 

(13) go(x) = { 

{ fidx Cdx 1 f?dx 2 ....f:«- 1 dx N f(x N ) , N>1 

for some a G R. Since g G C^°(R) we have F G «4.(R) which concludes the 
proof. □ 

The next theorem characterizes »4(R) as the set of Schwartz distributions 
with finite singular support which, in addition, are identical to a C°°(R) 
function on every open interval not intersecting its singular support. The 
theorem is not essential for the sequel and can be skipped. 

Theorem 3.2. A distribution F G £>'(R) belongs to A(M) iff: 

(a) The singular support of F is a finite set, and 

(b) If H C R is an open interval satisfying Hflsing supp F = 0, then there 
is a function g G C°°(R) such that the identity "F = g on H" is valid in the 
distributional sense. 



Proof. We will prove this result by using the definition of ,4.(R) that follows 
from Theorem 3.1. If F G T>'(M) satisfies the condition (a) above, we may 
set Vp =sing supp F which is then a finite set. Moreover, the restriction 
Fm\v f = h\v F for some / G C°°(R\V F ). We then have (F - f) R \ Vp = 
and so supp(F - /) C V F . Let A< F ) = F - f. Since is supported on 
the finite set Vf then (by a well-known textbook result [26]) A^' is a finite 
linear combination of Dirac deltas and their derivatives, i.e.: 
(14) 

N 

= £ A l F) wh ere A<£\x) =J2 C kHd^{x-w), N G N , C k (w) G C 

w£Vf k=0 

in agreement with eq.(8). Hence, F = / + A^) with A^) given by (8) 
and / G C°°(R\Vf)- To prove that F G A(R), by Theorem 3.1, we still 
have to prove that / G C^°(R). Let H =]xo,yo[ where xo,yo are two ar- 
bitrary consecutive elements of Vf- Then S n Vf = and from condi- 
tion (b) above F& = gs for some g G C°°(R). Since it is also true that 
F~ = fs it follows that fs = gs and so \im x ^ x + f^ = lvm x _^ x + and 

lirn - /" = lirn - g^ exist and are finite for all k G No- This proves 
that / G C£°(R) and so, by Theorem 3.1, that F G A(R). 
Conversely, if F G .A(R) then F is given by eqs.(7,8) for some finite set Vf 
and some / G C£°(R). Hence, sing supp F is a finite set. Moreover, on 
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an arbitrary open interval S C 1 satisfying 3Hsing supp F = 0, we have 
Fg = fw. Since / is infinitely smooth on 3 and f^ displays finite lateral 
limits at the boundaries of 3 (for all k € No), there is (by Whitney's ex- 
tension theorem [23]) a function g £ C°°(M) such that /= = gs and so also 
Fs = 55- This proves that F also satisfies condition (b) above. □ 

Before proceeding, let us make the following remark concerning the rela- 
tion between Vp and sing supp F: 

Remark 3.3. We will always assume that is a finite superset of sing 
supp F (but not necessarily identical to sing supp F). This is not contra- 
dictory with Theorems 3.1 and 3.2 (since sing supp F is always a finite set 
and we may always set A^ = for some w € Vp). The reason why we do 
not simply make Vf =sing supp F (although, most of the time, we assume 
this to be the case) is that the extra freedom will allow us to simplify the 
notation when deriving the explicit expression for the new product of dis- 
tributions. 

The aim now is to introduce a multiplicative product in the space »4(]R). 
We start by considering the (restriction to A(M) of the) product of distri- 
butions with non intersecting singular supports that was proposed by L. 
Hormander in [pag. 55, [12J]. 

Definition 3.4: The Hormander product • 

Let F,G G A(R) be such that V F fl V G = 0. Let {tl w cR,wGV F U V g } be a 
finite covering o/M by open sets satisfying w' £ Q, w , Vu) ^ w' £ VpU Vg- Let 
us also introduce the compact notation F w = Fq w , G w = Gq w to designate 
the restrictions of F and G to the set Q w . Then F ■ G is defined by: 



(15) F ■ G : (F ■ G)n w = F W G 



where F W G W denotes the usual product of a distribution by an infinitely 
smooth function. 

The Remarks 3.5 and 3.6 will show that the Definition 3.4 is consistent. 
They are included for completion. Remark 3.5 also provides the explicit 
form of F W G W . This result will be used in Theorem 3.7 to obtain an explicit 
expression for F ■ G. 

Remark 3.5. The product F W G W in eq.(15) is well defined. Indeed, let 
us write F,G € A(M) in the form given by eq.(7), i.e. F = / + and 
G = g + A^ G \ For each w £ VpUVo either F w or G w is a regular distribution 
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associated to a function in C°°(fi w ): 



(16) 



w e Vp 



F„ 



A (F) + f • 



/„GC»(^\H)n(7(o t 

Gui = G C°°(^u)) 



and 



w e Vg => 
where g w = gn u , and /„ 



Gi, 



A (G) + o ■ 



5w GC oo (^\{u;})nC7P(0 M ,) 

= fn w are the restrictions of g and / to the open 
set £l w . It follows that the product F W G W is well defined in the usual sense 
of a product of a distribution by an infinitely smooth function. We have: 



(17) 



F G 



g w A w '+g w f w , weV F 
w g Vg 



fw^w ~\~ gwfwi 



and so F W G W G V'(Q W ). 



Remark 3.6. The definition of F ■ G satisfies the prerequisites of Lemma 
2.5. Indeed, on an arbitrary open set S ww > = Q w nQ w > (with w ^ w') we have 
{F w G w )s ww , = {fgh ww , = (F w rG w ,) Sww ,. This is because E ww ,n{V F \JV G ) = 
0. Lemma 2.5 then implies that F ■ G is a Schwartz distribution, uniquely 
defined by eq.(15). 

We have not yet proved that F ■ G is independent of the particular cov- 
ering {£l w , w G Vp U Vq} used to define it. This result will be a corollary of 
the next theorem where we obtain an explicit expression for F ■ G: 



. Then 



Theorem 3.7. Let F,G G A(R) (such that V F H V G = $) be written in 
the form (7,8), i.e. F = f + J2 w£Vp A { W F) and G = g + Z we v G 
F ■ G is given by: 

F-G = fg+ 9 W ^L F) + E /- A - G) 

wGVf u>GVg 



(18) 



where g w and f w are C°°(R)- extensions of the restrictions g w = gn w and 
fw = fn w , respectively; and {ft w C K, w G Vp U Vg} is a finite covering of 
R satisfying the requisites of Definition 3.4- 

Remark 3.8. The extensions g w and /„, are C°°(IR)-functions such that, 
for each w G Vp we have g w (x) = g{x) if x G Q w ; and for each w G Vq, 
fw(x) = f(x) if x G fl w . That these functions exist follows from Whitney's 
extension theorem [241 [TO] and the fact that f,g£ C^ (M) and so, for each 
w G Vp (respectively w G Vg), the restrictions g w G C°°(0^) -cf. eq.(16)- 
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(respectively f w G C°°(Q W )) and all its derivatives display finite limits at 
the boundaries of Q w . Hence, the right hand side of eq.(18) is a well defined 
Schwartz distribution. Let us denote it by J. 

Remark 3.9. We now prove that J is independent of the particular fi- 
nite covering of M (satisfying the conditions in Definition 3.4) and of the 
particular extensions of g w and f w that were used to define it. Indeed, if 
{Q w , w G VfDVg} and {Q' w w G VfUVg} are two such coverings and g w , f w 
and g' w , f' w the associated functions then w ^ supp(/4 — f w ), G Vg and 
w <£ supp(<4 - g w ), \/w G V F . Hence: 

supp (g' w -g w )n supp A^f ) =0, Vw G V F 

and 

supp {f' w -f w )n supp A<£) = 0, y we V G 

and so 



(19) y, ^i F) + E /- A ^ G) 



/ j i>w W ' / j J w w 
w&Vf wEVg 







Therefore J is a well defined distribution that only depends of F and G. 

Proof of Theorem 3.7. It is straightforward to realize that the restriction 
of J to Q w is F W G W (given by eq.(17)). Since the restrictions to the open 
covering {£l w , w G Vf U Vg} define a unique global distribution (Lemma 
2.5) which, by definition, is F -G (cf. eq.(15)) the identity F -G = J holds. □ 

Two simple corollaries follow directly from the theorem: 

Corollary 3.10. F ■ G is independent of the particular open covering (Def- 
inition 3.4) used to define it. 

Corollary 3.11. If F,G G A(R) (such that V F nV G = 0) then F-G G A(R) 
and sing supp (F ■ G) C (V F U Vg). 

This concludes the study of Hormander's product. The next definition 
introduces the new product *, extending Hormander's product to the case 
of distributions with intersecting singular supports. 

Definition 3.12: The product * 

Let F,G G «4(K) . The multiplicative product * in A(M.) is defined by: 
(20) F*G= lim F ■ G e 

where G £ (x) = G(x + e) is the translation of G by e and the limit is taken 
in the weak sense. 
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The two following Remarks discuss the consistency of Definition 3.12. 
Then, in Theorem 3.15 we will prove that * is an internal operation in A(M.) 
and determine the explicit form of F * G for general F,G£ A(M). 

Remark 3.13. Since V F and Vq are finite sets of real numbers there 
is always a real number a > such that Vf D Vg^ = 0, Ve g]0, <t[; no- 
tice that Vg^ = {x — e, x G Vq] and so we may set, for instance, a = 
min ({\y — x\, x G Vf, y G Vg}\{0}) if this minimum exists and a = 1 if this 
is not the case (because Vf = Vq and $V F = 1)- Hence the product F ■ G e 
is well defined for e g]0, <j[. 

Notice that, in general, F G is not well defined since sing supp Ff) sing supp 
G 7^ 0. However, when that is the case, we will easily find (from Theorem 
3.15) that F -kG = F ■ G. 

Remark 3.14. The functional F-kG acts as: 

(21) <F*G,t>= lim <F-G e ,t> 

e— +0+ 

and if its domain is T>(M) - i.e. if the limit on the right hand side exists for 
all t € P(R) - then F * G is a Schwartz distribution (cf. Lemma 2.6). The 
next theorem will prove that this is the case for general F,Ge A(M) and 
furthermore that F*G G A(M). 

Before we proceed we need to introduce some notation. Let F,Ge A(M.) 
and (for notational convenience) let us set Vf = Vq =sing supp F U sing 
supp G. Let also N = U(Vjr) and Xk G Vf, k = 1, ..,iV, x^+i > be the 
array of elements of Vp (or Vq)- Hence, Vf = Vq = {xi, xn} and the 
distributions F and G can then be written in the form (cf. eqs(7,8)): 



(22) 



/ + ELiA (F) 



■j-'k 



where we set A^P = for Xk G Vising supp F and A^ = for Xk G 
y G \sing supp G. Moreover, f,ge C™ (R) n C°° (R\V F ) and so (cf. Remark 
2.3) they can be written in terms of a set of C°°(lR)-functions fk,9k', k = 
0,..,N: 

N-l 

f(x) = H(x 1 -x)f (x) + y^H(x-x k )H(x k+1 -x)f k (x) + H(x-x N )f N 



(23) 



k=i 



N-l 



g(x) = H(xi - x)g (x) + ^ H{x - x k )H(x k+1 - x)g k {x) + H(x - x N )g N 

k=l 

We can now state the new theorem: 
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Theorem 3.15. Let F,G € A(R). Then F * G e .A(R) and its explicit 
expression is given by: 

N 

(24) F * G = fg + £ [dk^S + A-i A i? 

fc=i 

for F, G written in the form (22,23). 

Proof. The first step is to obtain an explicit expression for F ■ G e (cf. 
Definition 3.12). Let then F,G be written in the form (22,23). It follows 
from eq.(22) that G e (x) = G(x + e) is given by: 

N 

(25) & = <f + £ A ^-e 

fc=l 

where g e (x) = g(x + e) and A^_ e (x) = A^(x + e) is a singular distribution 
with support on {xk — e}. We also have Vf = Vq = {xk, k = 1, .., N} and so 
Vg 6 = {xk — e, k = 1, .., A^}. Since there is a er > such that Vf H Vg« = 
for < e < cr (cf. Remark 3.13), we can use eq.(18) to calculate F ■ G e 
explicitly for all e <G]0, a[: 

(26) F ■ G e = fg e +Y, (5 £ ).A(f ) + Y /- A - G) 

W&Vf W&Vae 

N N 

fc=l fc=l 

where (<7 e )a; fe and /( Xfc -e) are the C°°(R)-extensions of the restrictions (g e ) Xk = 
(g e )n x and f(x k -e) = /n (xfc _ e )j respectively. As usual the open sets fi Xfc and 
ft/ Xk _ e \ belong to a finite covering {Q w c I; w G Vj? U Vg«} of R, satisfying 
the conditions stated in Definition 3.4. That is: i) U^ =l (Q Xk U Q/ Xk _ e -\) = R 
and ii) «/ ^ n w , Vw; w' £ (Vf U Vg«)- There are, of course, many possible 
coverings of this sort but the expression (26) is independent of the particular 
one used to define it (cf. Remark 3.9). One possibility is (let xq = — oo and 
xn+1 = +oo): 

(27) n Xk =]x k - e,x k+1 - e[ and Q Xk - e =]x k -i,x k [, ,k = l,..,N 
It then follows from eq.(23) that: 

(28) (g e h Xk =(g € k h Xk and (/)n (Xfc _ £) = (fk-i)n (xk _ e) , fc = l,..,JV 

where g\ is defined by = g k (x + e) and g^, were introduced in 

eq.(23). Since g k and belong to C°°(R) for all k = 0, .., N we may set (cf. 
eq.(26)): 

{9 e )n Xk =9k and (/)n (lSfc _ e) = > & = 
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Substituting this in eq.(26) we finally get an explicit expression for F ■ G € : 

N 

(29) F ■ = fgf + £ [s£ Ag) + A-i Agl, 

fc=i 

To obtain an explicit expression for F*G we still need to calculate the limit 
in eq.(21). Using eq.(29) we get for every t G V(R): 

( N N 

lim + <F.&,t> = £ lim </'9 £ 1 t>+^<#?,i>+^<A-iA£ 1 

I k=l k=l 

N 

> 



= lim + J f(x)g(x + e)t(x)dx + ^\im + <Ai F k \tg e k 

N 

(so) +E e l™ + < A £ ) ' re 4"-i> 

k=l 

N 

= / /(xMx)<(x)dx+x;<A(j),^ fc >+x;<Ag ! ),t/ fc _ 

J k=l k=l 

= <^ + E[^A£) + / fc _ 1 Ag)],t> 
fc=i 



where we defined t e (x) = t(x + e), /^(x) = fk(x + e) and used the fact 
that t,gk, fk £ C°°(R) and so they are uniformly continuous on an arbitrary 
compact interval. We conclude from eqs. (21,30) that F*G6 P'(R) and is 
indeed given by eq.(24). Moreover, (24) is obviously of the form (7,8) and 
so F*G G A(M). n 

We proceed by studying some of the properties of the new product: 

Theorem 3.16. The -k-product is (i) distributive, (ii) associative but (Hi) 
non- commutative, (iv) it reduces to the standard product of an infinitely 
smooth function by a distribution if either F or G belong to C°°(R) and (v) 
it reproduces the standard product of continuous functions. 

Proof. 

(i) Let us prove that the product is left distributive. Let F,G,J <G A(M). 
We have: 

(31) (F + G)* J= lim (F + G) ■ J e 

We now prove that if K G A(M) and Vf+g = Vf U Vg is sucn that Vf+g H 
V K = % then: 

(32) (F + G)-K = F-K + G-K 
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This identity is valid globally if it is valid locally. Let {£l w , w £ Vf+g U Vk"} 
be an open covering of M satisfying the conditions of Definition 3.4 for the 
• product of F + G by K. Then eq.(32) is equivalent to: 

(33) [(F + G) ■ K]a w = [F ■ K + G ■ K]q w , VweV F+G UV K 
which, in turn, follows from: 

(34) [(F + G)-K]n w = (F + G) W K W = F W K W + G W K W 

= [F.K]n w + [G-K]n w = [F-K + G-K]n w 

where we used the definition of the • product (eq.(15)) and the fact that 
for each w either F W ,G W e C°°(n w ) or K w € C°°(n w ). Making K = J e in 
eq.(32) and substituting in eq.(31), we get: 

(35) (F + G)*J= lim F • J e + lim G-J e = F-kJ + G-kJ 

which proves that the product is left distributive. Equivalently, one may 
prove that it is also right distributive, 
(ii) Let F,G,J € A(R). Then: 



(36) F = f+ A i F) > G = g+Y.^ ] and ^ = i+E A - 



w£Vp w€zVg w&Vj 



where f,g,j 6 C£°(M). To keep the notation simple let us redefine Vf = 
Vg = Vj =sing supp F U sing supp G U sing supp J. Let TV = JJVp and let 
Vf = {a?i, xn} where Xi < xt for i < k. Then F, G can be written as in 
eq.(22) with f,g given by (23) and: 



(37) J = j + Y, A 



N 

(J) 

k=l 



(F) 

(notice that Ax k = if x k G Vi?\sing suppF and equivalently for the 
other distributions) where j can be written in terms of C°°(lR)-functions 
j k ; k = 0, ..,N, as: 
(38) 



N-l 



j(x) = H(x 1 -x)j (x)+ ^2 H(x-x k )H(x k+ i -x)j k (x) + H(x - x N )j N (x) 



k=l 



A trivial calculation using eq.(24) shows that: 

N 



(39) {F*G)*J 

N 



-k J 



/^EK ) +/^») 
k=i 

fgj + £ {adk^ + fk-ijk^ff + Sfc-i/fc-iAi-j) 



fe=i 
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and also: 
(40) 



F*(G*J) = F* 



9J + +9*-i*g) 



k=i 



N 

= J 93 + £ (gkjktig + fk-dkAfJ + g k ^f k ^A^) 



k=l 



Hence, the ^-product is associative. 

(iii) Let F, G be given by eq.(22,23). We have from eq.(24) that: 



N 

(41) F*G-G*F = Y J [{9k ~ 9k-i) Ag? + {fk-i - fk) A<g> 

k=i 

and thus the ^-product is non-commutative. Notice nevertheless, that supp 
(F*G — G*F) CVpU Vq, is a finite set. 

(iv) It follows directly from eq.(24) that if F = f (or G = f) for some 
/ G C°°(R) then F -kG = fG (respectively, F*G = fF). 

(v) If F,G € (>A(M) n C°(R)) then, in the notation of eq.(22), F = f and 
G = g. Hence, from eq.(24) F * G = fg thus concluding the proof, n 

Hence A(M) is an associative algebra. To proceed we may endow A(M.) 
with the following bracket structure: 

(42) [F, G] = F ★ G - G * F 

the explicit form of the bracket being given by eq.(41). 



Theorem 3.17. The bracket (42) is a Lie bracket. 



Proof. The antisymmetric property of the bracket follows by construc- 
tion while the linearity and the Jacobi identity are inherited from the left 
and right distributive and associative properties of the ^-product, respec- 
tively. It is also easy to check that the bracket satisfies the Leibniz rule with 
respect to the *-product.n 

Finally we prove that: 

Theorem 3.18. Let 4- be the usual distributional derivative in T>'(R). 

The restriction ^ : A(M.) — > A(M) satisfies the Leibniz rule with respect to 
the -k-product. 

Proof. Let F, G G A(R). Then F*G G A(R) and its derivative acts 

as (F' = £F): 

(43) 

< (F*G)',t>= - < F*G, t' >= - lim < F-G\t' >= lim < (F-G € )',t> 
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Let us then prove that the product • satisfies the Leibniz rule. Let F,j£ 
A(M) be such that VfC\Vj = 0. Let {Q w , w € VfU Vj} be an open covering 
of R satisfying the conditions of Definition 3.4. Locally we have: 

<=► K F -JhJ=[F'-J] nw + [F-J'] nw 
^ '* {F w J w ^j — F W J W -\- F W J W 

which is true because either F w or J w belongs to C°°(Q W ). Hence the • 
product satisfies the Leibniz rule locally and thus also globally. 
Substituting this result in eq.(43) we get: 

(45) <(F*G)',t> = lim < F ■ (G e Y + F' ■ G e ,t > 

= lim < F ■ (G'Y , t > + lim < F' ■ G\ t > 
= < F*G',t> + < F'*G,t> 

where we took into account that the translation and the derivative operators 
commute in V(M). Hence, (F * G)' = F *G' + F' *G, which concludes the 
proof. □ 

Corollary 3.19. The algebra (*4(R), +, *) is an associative (but non- 
commutative) differential algebra of generalized functions. 

This concludes the study of the properties of the ^-product and of the 
associated algebra ^4(R). As a simple example, let us calculate the derivative 
of H(x) *H(x). We have: 

— (H(x) * H(x)) = 6(x) ★ H(x) + H(x) ★ S(x) = S(x) * H(x) = 6(x). 
ax 

which of course, is consistent with the fact that H(x)-kH{x) = H(x). Notice 
that multiple products of the Heaviside step function are sometimes used 
to exemplify a simpler version of the Schwartz impossibility result [3]. One 
typically finds that differentiating the equation H n (x) = H(x) for different 
values of n leads to contradictory results. This is not the case if H n (x) is 
calculated using the product *. The key point, as one can easily check, is 
that the * product is non-commutative. 

To finish this section, let us point out that there are other associative but 
non-commutative products, related to the one introduced in eq.(20), and 
also yielding a Lie bracket structure in »4.(R). For instance: 

F * 2 G = lim F ■ G~ e 

(46) F* 3 G = lim F e ■ G 

F* A G ee lim F~ e ■ G 
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which are all related to the original product introduced in (20): 

(47) F* 2 G = F* 3 G = G* A F = G*F 

Finally, we can also define a commutative product (or symmetric bracket) 
through the prescription: 

(48) F* 5 G = ^(F*G + G*F) 
which however is not associative. 



4. Linear differential equations with distributional 

coefficients 

In this section we consider the second problem presented in the Introduc- 
tion. Let: 

n 

(49) J>V> (<) =/, 

i=0 

be an ordinary, linear differential equation defined on R. The coefficients 
are the complex valued functions a. L G C 00 (M) and a n {x) ^ 0, Vx € R, 
and the non- homogeneous term is / G C°°(M). As usual, '0 denotes the 
i-th order derivative of tp. Let ipu be the general solution of eq.(49). The 
subscript U stands for unconfined (meaning the ipu is the solution of (49) 
on the entire real line). The generalized solution of eq.(49) (i.e. the solution 
in Z>'(R)) coincides with its classical solution (i.e. the one in C°°(R)) |14| . 
The designation tpjj will be used to denote both the distribution in 2?'(R) 
as well as the associated infinitely smooth function. 

Let us also consider an object of the form ipc = H(x)ipjj "confined" to the 
positive axis and reproducing the "unconfined" general solution for x > 0. 
Such a function is not a generalized (nor classical) solution of eq.(49). The 
aim of this section is to derive a new differential equation that displays the 
generalized solution (and only the generalized solution) ipc- 

Notice that we are focusing on the particular case of Problem 2 (stated in 
the Introduction) where $7 = R + . This is to keep the formulation as simple 
as possible. The extension to the general case where £1 =}a,b[, a < b G 
R U { — do, +00} is straightforward. 
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Let us then proceed. As a first step we substitute ipc m eq.(49) and find 
the correction terms: 



(50) £ofl$ = ^aAj2Q)H^4' j) 



i=0 i=0 [ j=0 

n n i 

= ^ H 4 +£«*Eooff (i) Wr ) 

2=0 1=1 j = l 

i>j=i 

Remark 4.1. Equation (50) is not a closed equation for ipc, because it 
displays a distributional non-homogenous term which is dependent on the 
particular solution of eq.(49). For each solution of eq.(49) (associated to a 

complete set of boundary conditions ip^\xo),i = 0, ...,n — 1 given at some 
xq £ M + ), this term can be re- written as a combination of Dirac deltas and 
their derivatives: 

n n i— 1 

£ {))a t {x)H^{x)4~ j \x) = J2J2< x ^ { u(°) d(i ~ 1 ~ k) ( x ) 

i>j=l i=l fc=0 

and upon substitution in (50) we obtain a linear and ordinary differential 
equation, with a distributional non-homogenous term, that for a particular 
set of boundary conditions, displays a confined solution. However, the gen- 
eral solution of this equation is still not confined to M + , but exhibits instead 
a jump discontinuity at x = [H]. More details about ordinary differen- 
tial equations with distributional non-homogeneous terms can be found in 
[Hi [TBI [23] . The related issue of distributional solutions of ordinary differen- 
tial equations is studied in [TTJ, [151 [2S] 

Our task here will be to re- write eq.(50) exclusively in terms of ipc so that 
its general solution (and not only a particular solution) is of the form Hipjj. 
With this in mind, we prove the following theorem: 



Theorem 4.2. // ipu € C°°(R) and ipc G A(M) is such that ip c = 4>u 
on M + then, in the distributional sense: 

(51) £ ^a iH ^4~^ = £ ())a lH ^ *4~i\ 

i>j=l *>i=i 

Proof. The most general distribution ipc £ A(M) satisfying ipc = ipu on 
M + is given by: 



(52) 



ipc = Hi^u + F 
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where F E A(M) is such that supp F C R U {0}. Hence we have (for all 
i>j = 

<H^K^" j \t> = <H {j) *[H^u + Ff~ j \t> 

= < *[H* fa] ,t> + < * irfi-i) ) t > 
= < ★ ^ ( r~ j) , t > + hm < ■ (F^-V , t > 

e^0+ 

(53) = <H®il>§- j \t> 



where in the last step we took into account that for j > 1, the sets supp 
HW = {0} and supp (F e p z ~^ C] — oo,— e] are disjoint. The result (53) is 
valid for all £ > j = 1, ..,n and so eq.(51) holds. □ 

We can now re- write eq.(50) exclusive in terms of tpc- 



n n 

(54) £ = If / + £ QH^) * ^ 

i=0 i>j=l 

which seems to be a suitable candidate for the differential equation we are 
looking for. That ipc = Htpjj is a solution of eq.(54) is obvious by construc- 
tion. However, we still have to show that is its only solution: 

Theorem 4.3. Let ipjj be the general solution of eq.(49). Then, the general 
solution of eq.(54) in A(M) is given by tpc = Htpu- 

Proof. For a generic ipc € A(M) the support of the second term on the 
right hand side of eq.(54) is {0}. Hence, on R + and R~ , eq.(54) reduces to 
eq.(49) and to the homogeneous equation associated to eq.(49), respectively. 
We conclude that any solution of eq.(54) will be of the form: 



(55) Vc = A + + 



where supp A C {0} (i.e A is a linear combination of the Dirac measure and 
its derivatives), H-(x) = H(—x) is the reversed Heaviside step function H 
and -0_ , tp + £ C°° (R) are the distributional (which coincide with the classical 
[14] ) solutions of Y^=Q a i' l l ) = an d °f ec L-(49), respectively. Substituting 



A MULTIPLICATIVE PRODUCT OF DISTRIBUTIONS AND APPLICATIONS 21 

(55) in eq.(54) we get: 

i=0 i=0 j=0 

= Hf+ Q)aiH^4- j) 
i>j=i 

n n 

J>A« + J> (H^{x)+H^) 

i=0 i=0 

n n 



Separating the terms that involve the delta distribution or its derivatives 
from those that do not, we get: 

n 

(57) Yl a * { H -^ + ^+) = F ^ 

eiu^o^ 0*0 = /(*)> *>° 

, E?=o a i( x W-( x ) = , x<0 
(confirming what we already knew about and: 

n n 

(58) J> t AW - £ (jOoifl-^tf - ^ = 

i=0 *>j=l 

In eq.(58) the term proportional to the highest order derivative of 5 is given 
by a n A^ (if not zero, a n A^ is proportional to at least 5^). This term 
cannot be cancelled by a combination of lower order derivatives of 5. Since 
a n ^ we get A*™) = and so A = 0. Hence, eq.(58) reduces to: 

n n n 

(59) £ {)) ai H^:- j) = o^Yl EG-)**- = 

i>jr' = l j = l j=j 

From eq.(57) we find that Y^i=o a i^-' l l J<y - = 0, Vx G 1R and by adding this 

last equation to eq.(59) we finally get: 

(60) 

n n n w n 

Y^y^.ir 1 c' 1 = o ^ ^> J- (h^_) = o J>*« = o 

j'=0 i=j i=0 i=0 



where (f> = H_\[)_ satisfies (by construction) the boundary conditions 4>^ 1 \zq) ■ 
0, i = 0, ...n — 1 at an arbitrary zq > 0. For these boundary conditions the 
unique solution is (f> = 0, which implies that ijj-(x < 0) = 0. 
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Assembling all these results, we finally get: Vc = Hvp + , where ^M- satis- 
fies eq.(57) (or equivalently eq.(49)), which proves the theorem. □ 

We conclude that eq.(54) provides a global formulation for the differential 
problem described by eq.(49) and confined to the positive axis, i.e. it yields 
a solution of Problem 2 for the case where f2 = R + . The extension to the 
more general case where fi is an arbitrary interval is straightforward. 
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